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Abstract. We consider upper bounds for the approximation er- 
ror M\g(X) — g(X)\ p , where X and X arc random variables such 
that X is an approximation of X in the L p -norm, and the function 
g belongs to certain function classes, which contain e.g. functions 
of bounded variation. We apply the results to the approximations 
of a solution of a stochastic differential equation at time T by the 
Euler and Milstein schemes. For the Euler scheme we provide also 
a lower bound. 



I. Introduction 

Convergence schemes for the solutions of SDEs are rather well known. 
Let X be the solution of the one- dimensional equation 

dX t = a(t, X t ) dW t + b(t, X t ) dt, X = x , 

where W is a standard Brownian motion, t G [0, T], and a and b satisfy 
certain assumptions specified in Section 2. P.E. Kloeden and E. Platen 
[12] have showed that any order of strong convergence can be achieved 
by the strong Ito- Taylor approximations, i.e. for any order 7 > there 
exists a scheme X w corresponding to a partition tt of the interval [0, T] 
with mesh size \tt\ such that 



E sup \X t -X?\ < C|tt| 7 . 

\0<t<T / 

The simplest and most common examples are the Euler scheme X E and 
the Milstein scheme X M , which have the order of strong convergence 
0.5 and 1, respectively. 

Lately this topic has been considered by N. Hofmann, T. Miiller- 
Gronbach and K. Ritter in [7, 8], Hofmann and Miiller-Gronbach in 
[9], and Miiller-Gronbach in [14, 15]. In these papers they cover errors 
with respect to both global and pointwise error criterions. The latest 
result concerning the pointwise error is due to Miiller-Gronbach [15], 
where the author defines certain classes of convergence schemes and 
finds optimal (adaptive) schemes for each class. 
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The pointwise weak error E/(Xy) — E/(X£) is also considered by 
several authors, e.g. Kloeden and Platen [12], under certain smooth- 
ness conditions on /. In the case of the Euler scheme the smoothness 
conditions were relaxed by V. Bally and D. Talay in [2], where / was 
required to be only measurable and bounded. This was done under 
a Hormander type condition for the infinitesimal generator of X. A 
recent contribution to this field is given by Moon et al. in [16]. 

This paper considers the error of the type E l \g(X T ) — g(X^)\ p . The 
question is motivated by discretization schemes for BSDEs. The termi- 
nal condition q(Xt) is approximated by g(X?f), and L p -estimates for 
the difference g{Xx) — g{X?f) are required. If g is Lipschitz, this error 
returns immediately to the error of the underlying scheme. Therefore 
the aim of this paper is to give results for relevant non-Lipschitz func- 
tions. It is proved that if approximations (X^) t e[o,T\ satisfy 

\\x T - xR p < c\ kr 

for some 7 > and all 1 < p < 00, then 

E\g(X T )-g(X^<C 2 p \irr £ 

for any < e < 7 and for any g in a special class of functions. This 
class contains functions of bounded variation, polynomials and jump 
functions controlled by the tail distributions of X T and and there- 
fore by Minkowski's inequality the result is true for any sum of these 
three types of functions. 

Finally, for the Euler scheme a lower bound is presented indicating 
that the error under consideration can not converge faster than |7r| 1//2 . 
This is done by giving an example, namely the geometric Brownian mo- 
tion, for which the lower bound is obtained. There is still a difference, 
although arbitrarily small, between the upper and lower bounds, and it 
remains open whether the rates could be equalized. These results are 
achieved under certain conditions on the SDE, including the existence 
of a bounded density for the solution Xt- 

The structure of the paper is such that the assumptions that hold 
throughout the paper are collected in Section 2. Sections 3-5 contain 
the main results concerning upper bounds. The first result is given in 
Section 3, where Theorem 3.4 clarifies the convergence rate for indica- 
tor functions. This is then applied to the Euler and Milstein schemes 
in Theorem 3.7. The result is then extended to functions of bounded 
variation in Theorem 4.3 in Section 4, and applied to the Euler and 
Milstein schemes in Theorem 4.5. Another extension is developed in 
Section 5, where the result for the function class Q v ^ is given in The- 
orem 5.7. The class G P .<p is then analyzed in Section 6, including the 
result that it contains all polynomials in Theorem 6.3. An application 
to the Euler and Milstein schemes is presented in Corollary 6.6. Section 
7 contains a lower bound for the convergence, stated in Theorem 7.2. 
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Finally, a proof of a Theorem from the book of Bouleau and Lepingle 
[1] is presented with explicit constants in Appendix A. 

2. Assumptions 

We fix a terminal time T > and suppose that (W t )te[o,T] is a stan- 
dard one-dimensional Brownian motion defined on a complete filtered 
probability space (Q, J 7 , P, (J 7 t)te[o,T\) , where the filtration is the aug- 
mentation of the natural filtration of W and T = Tt- 

We consider a diffusion process X, which is a solution to 



(2.1) 



dX t = a(t, X t ) dW t + b(t, X t ) dt, 

Xq = Xq 



with xq G R and a, b : [0, T] x R — > R. We assume that a, b G 
C 0,1 ([0,T] x R) and for / G {c, 6} there exist constants Ct and a > \ 
such that 

(i) |/(t,x)|<L7 T (l + |x|), 

(ii) \f(t,x)-f(t,y)\<C T \x-y\, 

(iii) |/(t,x)-/(s,x)| <C T (l + |x|)|t-s| a . 

Assumptions (i)-(iii) imply the existence of a unique adapted strong 
solution X of the SDE (2.1), see e.g. [11, p. 289]. For estimates 
concerning the Milstein scheme we assume in addition that 

(iv) the state derivatives || and || satisfy the condition (ii). 
Moreover, we assume that 

(v) X T has a bounded density. 

Remark 2.1. Assumption (v) is satisfied (see [6, p. 263]), if we as- 
sume that a, b G C£°([0, T] x R) and er satisfies the uniform ellipticity 
condition, i.e. there exists a constant f3 such that 

a(t,x) > (3 > for all (t,x) G [0,T] x R. 

Another sufficient condition is given by Caballero et al. in [5, Theorem 
2]. They assume that a and b are C 2 in x, the second derivatives 



have polynomial growth, the functions \a(0,x) 
\b x (t,x)\ are bounded, and 

-P()/2 N 



\a x (t,x)\, 1 6(0, x) I and 



E 



/ a(s,X s ) 2 ds 
Jo 



< oo 



for some p > 2 and for all t G (0, T]. Then there exists a continuous 
density /x t of X t such that for all p > 1 

'1/2 



/x t (x)<C7 p 



a(s, X s ) ds 



for some constant C p > 0. 
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Denote by n a partition = t < ti < . . . < t n = T of the interval 
[0,T], and let 

\n\ = max \t i+1 — tA 

0<i<n 

be the mesh size of ir. Moreover, denote an approximation of X cor- 
responding to 7r by X n . Two such approximations are the well known 
Euler and Milstein schemes. 

Definition 2.2 (Euler scheme). Let X E be the Euler scheme relative 
to n, i.e. Xq = Xq, and for i — 0, . . . , n — 1, 

X t E i+i = X* + b(U, X t E J(t i+1 - U) + a(U, X t E )(W ti+1 -W ti ). 

Given the values at the partition points, we also define the Euler scheme 
in continuous time by setting 

Xf = Xf k + a(t k ,Xf k )(W t - W tk ) + b(t k ,X t E k )(t - t k ) 

for t G (tk,tk+i)- This can be written using the integral representation 

+ / J2 b ^ k ' X 0x[t k ,t k+1 ){u)du 

Jo k=0 

for t G [0, T] a.s., which coincides a.s. with the discrete scheme at the 
partition points. In this paper we need the continuous time version 
only for the Euler scheme. 

Definition 2.3 (Milstein scheme). Let X M be the Milstein scheme 
relative to n, i.e. X^ = x , and and for i — 0, . . . , n — 1, 

x t M +i =X% + b(t t , X t M )(t l+1 - u) + <t(U, xff)(w ti+1 - W u ) 

+ \a{U,X^)a' x {U,Xt;){{W ti+1 - W u f - (t i+1 -t,)). 

We remind that the Euler and Milstein approximations of X^, de- 
noted X E and X^ , always depend on the corresponding partition it. 
This is omitted from the notation for simplicity. 

These assumptions hold throughout the paper. 

3. Indicator Functions 

3.1. General Approximation. Suppose that we have a probability 
space (f2, J 7 , P) and two random variables X, X : Q — > R. Consider X 
to be an approximation of X in the L p -norm. First we find an estimate 
for the quantity E|g(X) — g(X)\, where g is an indicator function. 

Definition 3.1. Recall the non- increasing rearrangement of a random 
variable X, defined by X* : [0, 1] — > R U {+oo, -co}, 

X*(s) := inf{c G R : P(X > c) < s}. 
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Here we use the convention that inf = oo. 

Remark 3.2. Definition 3.1 is slightly different from the standard non- 
increasing rearrangement as defined e.g. in [4], where the absolute 
value of the function X is taken. However, by analoguous arguments 
we can show the following properties: 

(i) X*(l) = — oo, X*(0) = oo if X is not essentially bounded and 
X*{s) G R for s G (0,1), 

(ii) X* is right-continuous, 

(iii) X* has the same distribution as X with respect to the Lebesgue 
measure on [0,1], 

Definition 3.3. Denote the minimal slope of the function X* from the 
level K by d x ■ R -> [0,oo), 

se[o,i] [ \s — a(K)\ J 

s^=a(K) 

where 

a(K) = P(X > K). 

Theorem 3.4. Suppose that X is a random variable. Then the follow- 
ing assertions hold: 

(i) If X has a bounded density f x , then for all K G R, all random 
variables X and all < p < oo we have 

n X [K,oo)(X) - X [K,oo)(X)\ < 3D X (K)^ 
where 

Dx{K):= dxjK) e ^ SUpfx] - 
Moreover, the power of the L p -norm is optimal, i.e. if 

(3.1) E\ X[ k,oo)(X) - X [K,oo)(X)\ < C(X,K,p) 

for all random variables X, then the power can not be re- 
placed by a power q such that < q < oo. 

(ii) // there exists po > such that the formula (3.1) holds for all 
Po < P < oo, all K G R and all random variables X , and 
there exists B x > such that C(X,K,p) < B x , then X has a 
bounded density. 

Proof. Let us first show (i). Fix K G R and < p < oo, and let X be 
a random variable such that 

E|x[x,oo)P0 -X[K,oo){X) \ =e 

for some e G (0, 1]. Define e 1 := P(X > K, X < K) and e 2 := P(X < 
K, X > K), so that e = E\ + e 2 . Denote by a the number a(K) 



X -X 



X -X 
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introduced in Definition 3.3 and notice that a — S\ > and a + e 2 < 1. 
Now 

E|X - X| P > E|X - X\ P X{X>K,X<K}\J{X<K,X>K} 

— E|X — -ft'| P X{x>_ft:,x<E'}u{x< J ft',x> J K'} 

= E|X - + E|X - Xr X{ x<^>^ } - 

Since X has a bounded density, we can find a number Co G [X, oo] such 
that P(K < X < c ) = ei, thus also \{K < X* < c }\ = e x . Note 
that c may not be unique. But {K < X < c } is a set of probability 
ei where E|X — K\ p xa is minimized over all A C {X > .fT} with 
P(-A) = £i, which implies that 

E|X - X| P X {X >^ <K} > E|X - Xr X[x , Co) (X) 
= / |X*( S )-Xr X[XiCo) (X*(s))rf S = f |X*( S )-Xp>d S 

</[0,l] Ja-ex 

~ Jo P+1 
and by similar arguments 

" a+£2 d x (KYe p 2 +1 



p+l 



i>a+E2 

E\X - K\v X{x<K , x > K} > / \X*(s) -K\*ds> 

J a 

Thus 

(3.2) m -x r >M!DM^n > Mp^. 

K J 1 1 ~ p+l ~ 2P(p+l) 

Now the equation (3.2) gives 

nX[K,oo)(X) - X [K,oo)(X)\ 

P 

< 2 *o> + i)* (_L_p( E |x-xr)^ 

By elementary computations we can show that 

2?+r(p+ i)irr < 2 e A < 3, 
and keeping in mind the definition of D x we can write 



E| X [*,oo)P0 - X^,oo)(X)| < 3DxW^ 



X -X 



V 

p+l 



Using the definition of X* and the boundedness assumption for the 
density of X we see that l/d x (K) < sup fx- 

Moreover, the power -^y of X — X is sharp. To see this, we 

p p 

construct an example where the lower bound given by equation (3.2) 
is achieved. Suppose that Q = [0, 1] is equipped with the Lebesgue 
measure, K — \ and e < 1. If we take X(uo) = u, then X has a 
bounded density and d x {\) = 1. Now define 
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X 



x. if ^ G [o, i -§)u + 

X + §, if wG [± 
X- 



e 
2' 



2 2' 
L2 2' 2J' 

if we a, i + §]. 



Then 



E|X-X| p 



E 



X[i- f ,i +f] P0 



2p 



2p 



which coincides with the lower bound in equation (3.2) up to the con- 
stant. Hence the power of 



X -X 



can not be increased in the 



assertion (i). 

Now we verify (ii). Let 5 > and choose X = X — 5. Then 

E|X[K,oo)P0 - X[X,oo)(^)| 

= P(X >K, X -5 <K) + P(X <K, X -5>K) 
= P(K <X <K + 5), 

so that by assumption we get, for p > p , that 

P(K < X < K + 5) < B x (E|5| p )^r < B X S^. 

We let p go to infinity and conclude that 

P(K <X <K + 5)< B x 5. 

Let N C R be a null set with respect to the Lebesgue measure and let 
e > 0. Since the Lebesgue outer measure of N is also zero, we find a 
sequence (Ij) of open intervals such that N C (J Ij and ^2 \Ij\ < e. Let 
Cx be the law of X. Then we have 



C x ((a,b)) < MM)) < S x |fe-a| 



and 



\ j / 3 j 

This implies that Cx(N) = 0, so £x is absolutely continuous with 
respect to the Lebesgue measure. By the Radon-Nikodym theorem 
there exists a measurable function / : R — > [0, oo) such that 

£ X (M)= f f(x)dx 

J M 

for all measurable M C R. Moreover, / is integrable since £x(R) = 1- 
Define a function $ : R — > [0, 1] such that 

$(t) = f f(x)dx. 
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Then by [18, Thm. 8.17] we have that $'(0 = /(*) a - e - in ^ 0n the 
other hand, we have that 

# . . , $(t + /i)-$(t) , „ 
&(t) = lim -i f ^ < lim — = B x a.e. in R, 

h^O h h^O h 

because <&(t + h) — = Cx((t,t + h)). Therefore we conclude that 
/(*) < B x a.e. in R. □ 

Remark 3.5. By considering complements of the intervals in the indi- 
cator functions and the random variables —X and —X, we have corre- 
sponding results for the functions X(k,oo), X(-oo,k] and X(-oo,k)- 

As an immediate consequence of Theorem 3.4, we can derive 

Corollary 3.6. Let X be the solution of the equation (2.1), K G R 
and < p < oo. Let X T have a bounded density and suppose that XJf 
is an approximation of X T such that 

\\X T — Xt\\ p < C p \n\ e 

for some 9 > and some constant C p > 0. Then for all K G R we 
have 

nX[K,oo)(X T )-X[K,oo)(XZ)\ < ZDxAK^cf 1 |tt|& . 

3.2. Euler and Milstein Schemes. Now we can apply the results of 
Section 3.1 to the Euler and Milstein schemes: 

Theorem 3.7. For any < e < 1/2 there exists a constant C £ > 
such that for all K G R we have that 

nX[K,oo)(X T ) - X[K,oc)(X*)\ < (D Xt (K) V y/D XT (K))C e 

and for any < e < 1 there exists a constant C' £ > such that for all 
K G R we have that 

nX[K,oo)(X T )-X[K,oo)(X™)\ < (D Xt (K) V VDx T (K))C' £ In] 1 - 6 . 

Proof. Let 1 < p < oo. Then for the Euler scheme we have by Theorem 
A.l in the Appendix that 

II E\\ - 

\\X T - X T || p < C p I 7t 1 2 , 
i.e. the assumption of Corollary 3.6 is satisfied with 6 — \. Thus 

(3.3) E\ X[ k,oo)(Xt) - X[K,oo)(X$)\ < ZDxAK^cf 1 H*T . 

Similarly for the Milstein scheme we have by [13, Proposition 1, p. 140] 
that 

\X T - X?\ p < C' p \ir\ , 
which gives the assumption of Corollary 3.6 with 9 = 1, and therefore 

(3.4) E\ X[ k,oo)(Xt) - X[K,oc)(X¥)\ < SD Xt (K)^(C' p )^ . 
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The claim follows in both cases by choosing p such that p — (6 — e)/e, 
where < e < 6, and noticing that for any a > we have ap+ ! < a V 
y/a. The constant 3 and the constants coming from the approximation 
schemes are included in C £ or C' £ , which now depend on e through the 
choice of p. □ 

Since we have information about the constant C p in Theorem A.l, 



o 2 

the Euler scheme: 



i.e. C p = e Mp , we can write an extended version of Theorem 3.7 for 



Theorem 3.8. Let K e R. Then there exists m e (0, 1) such that for 
\ir\ < m we have 



2+M 



E\x[K,oo)(X T )-X[K,oo)(X^\ < (D Xt (K)V^D Xt (K))\tt\~ 2 (-'oskl)^ , 

where the constant M = M(xo, T, Ct) G (0, oo) is taken from Theorem 
A.l. 

P 

Proof. By Corollary 3.6 and Theorem A.l, using 0/P+ 1 < a V y/a for 
a > and p > 1 , we get 

E|x[x,oo)(X T ) - X[x,oo)(X#)| < 3D Xt (K)^C^ |tt| W 



(3-5) < (,Dx T (^) V v / ^x T W)3e Mp 'FFr |tt| W 

< (£x T W V y/D XT (K))3e Mp2 \ir\ W 

for all p > 1. Now choose p such that 

4p(p + l) 2 = - log \tt\ 

for | -71" ] < to with to = e~ 16 .This gives p 3 < — log | tt | and p 2 < 
(— log |7r|) 2 / 3 . Thus we have 

e Mp 2 < e M(- log|7r|) 2 / 3 = | 7r |-M(-log|7r|)-V3 

and 



1 / p 



< J(-log|7r|)V3-i = (_l og | 7 r|)-V3. 



2(p + l) V _1 °gl 7r l 
Using these we get 

„ p 9 1 1 1 1+M 1 2+M 

3 e ^p 2 i^i^fft) = 3 e M P 2-27^+17 < 3 ( -io gW )i/3 < ^2 { - log]7rl) v3 ? 
where in the last step we used the inequality 



3 1 7T I (-logM) 1 ^ < I 

for I -7T I < to. Now we come back to equation (3.5) and conclude that 

n X [K,oo)(X T ) - X[ K,oo)(X E )\ 



1 2+M 

<(%(^)Vv^(^))H 5 (-iogM)i/3 



□ 
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4. Functions of Bounded Variation 

From Theorem 3.4 we deduce the same error for functions of bounded 
variation, up to a constant. Let us first recall the definitions of the 
spaces BV and NBV . 

Definition 4.1. Let 

N 

T f (x) := swp^2\f(xj) - f(xj-i)\, 

where the supremum is taken over N and all partitions — oo < xq < 
Xi < . . . < rrjv = x, be the total variation function of /. Then we say 
that / is a function of bounded variation, / e BV, if 

V(f) := lim T f (x) 
is finite, and call V(f) the (total) variation of /. 

Definition 4.2. Let NBV be the set of functions / € BV such that 
/ is left-continuous and f(x) — > as x — > — oo. 

4.1. General Approximation. 

Theorem 4.3. Suppose that X and X are random variables and X has 
a bounded density. If g e -BV and 1 < p < oo ; £/ien /or any 1 < g < oo 

E|o(X) - g(X)\P < 3 p+1 (sup f x )& V{g) 

Proof. First we show the result for functions g G NBV . By [18, Thm. 
8.14] there is a unique signed measure \i such that 

g(x) = fi((-oo,x)) and \/j,\((-oo,x)) = T g {x), 

where is the total variation measure of //. We consider the Jordan 
decomposition of /i , i.e. fj, — /ii — fj,2, where /ii = |(|//| + //) and 
A*2 = |(| A* | — are positive measures. Then |//| = //i + // 2 , and all 
three measures /ii and /i2 are finite since |/i|(H) = V(^) < oo. 
Thus we get 

g(x) = n((-oo,x)) = / X(-oo,x)(z) dfi(z) = / X(*,oo)(z) <^(*)- 

J R ./R 

Now by Theorem 3.4 and Remark 3.5 we get 



X -X 



9{X)-g{X) 



R 



X(z,oo)(X) dfi(z) - / X (z,oo)(X) dfi( 



R 



R 



X(*,oo)P0 - X(z,oo)(X) dfi(z) 
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< 



'R 



X(z,oo){X) - X(z,oo){X) d\n\(z) 

X{z,oo){X) - X(z,oo){X) d\fi\{z) 

< 3?(sup fx)^V(g) 



< 



R 



A- A 



P(9+1) 

9 



which completes the proof for functions in NBV . 

Next, let g be an arbitrary function in BV . By [18, Thm. 8.13], there 
exists a unique function g G NBV and a unique constant c£R such 
that g(x) = g(x) + c at all points of continuity of g, with V(g) < V(g). 
Also by [18] we know that g can have only countably many points of 
discontinuity, so define Uji^Oj} to be the set of these points and let 
Xj := g(a,j) — g(aj) — c. Then we can write 

g(x) = g(x) + c + A(x), 

where 

oo oo 

A (^) : = ^2*jX{ aj }(x) = J2 X i (X(-oo, aj ](x) - X(-oo, aj )(x)) ■ 

3=1 3=1 

We define a measure 

oo 

y = E A A> 

3=1 

where S a is the Dirac measure in a. Again by [18] we know that g(aj—) 
exists, so we have g(aj) + c = g(a,j—) and 

oo oo 

M(R) = \ x i\ = E \9^i) - 9(aj-)\ < V(g). 

3=1 3=1 

Now we can write 



A(x) 



R 



X(-oo,z](x) - X(-oo,z)(x) dv{z). 



and compute, similarly as in the NBV case, that 



ApT) - A(A) 



+ 



< 


In 


V 






-X(- 



\X(-oo,z](X) - X{-oo,z]( X )\ d W\( z ) 



R 



< 



R 



X(-oo,z] (A) - X(-oo,z] (A) d\u\ (z) 



+ 



R 



) (A)-X(-oo, z) (A) d\u\(z) 



< 2-3p(snpf x )^ T ^V(g) 



X-X 



p(«+i) 

9 
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This, combined with the NBV result, implies that 



9(X)-g(X) 



< 



~9{X)-~g{X) 



< 3 • 3p(supf x ) pi " +1) V(g) 
which gives the statement. 



g(X) - g(X) + A(X) - A(X) 
A(X) - A(X) 

X -X 



□ 



As in Corollary 3.6 for indicator functions, we can now write an 
analoguous statement for functions of bounded variation: 

Corollary 4.4. Let X be the solution of the equation (2.1), 1 < p < oo 
and g G BV . Suppose that X T has a bounded density, 1 < q < oo and 
X£ is an approximation of X T such that 

\X T - X^\ q < C q \n\ e 

for some 9 > and some constant C q > 0. Then 

E\g(X T ) - g(X})\* < 3 p+1 (sup/ XT )^ Vigfct 1 M^ 1 ■ 
4.2. Euler and Milstein Schemes. 

Theorem 4.5. Let g G BV and 1 < p < oo. Then we have for 
< e < 1/2 that 

E\g(X T ) - g(X$)\r < 3 p (su P / Xt V V^plZ)V(g) p C £ \n\^ £ 
and for < e < 1 that 
E\g(X T ) 



|l-£ 



g(X™)\ p < 3 p (sup/x T V ^/sup f XT )V( g yC' E \n\ 

where C e and C' £ depend on e and the constants of the corresponding 
schemes. 

Proof. The statement follows from Corollary 4.4 with arguments similar 
to the proof of Theorem 3.7. □ 

For the Euler scheme we can write an extended version corresponding 
to Theorem 3.8: 

Theorem 4.6. Let 1 < p < oo and g G BV . Then there exists m > 
such that for \ir\ < m we have 

2+M 



E\g(X T ) -g(X*)\ p < 3 p (sup/x T V y/^f^)V(g) p , 
where M is the constant in Theorem A.l. 

Proof. By Theorem A.l and Corollary 4.4 we get for 1 < q < oo that 
E\g(X T ) -g{X$)\* < 3 p+1 (su P / Xt V V^Pf^)V(g) p e Mq2 ^ W ' 
< 3 p+1 (sup/ XT V ^^hT^Vigye^ 2 



(9 + 1) 



2(9 + 1) 
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which by the arguments in Theorem 3.8 implies that 



1. 2+M 

E\g(X T ) -g(X*)\* < 3 p (su P / Xt V y/^f^)V(gY i-wW 



□ 



5. Extension 



Next we extend the result into a function class, to be called Q Ptip) that 
contains e.g. all polynomials. The main result is given in Theorem 5.7 
and the class Q V)ip is analyzed in Section 6. 

Definition 5.1 (Bump function). Let cp : R — > R be a function such 
that < if(z) < 1 for all z G R, ip is increasing in (— oo, 0] and 
decreasing in (0, oo), and 

lim ip(z) = 0. 

\z\— >oo 

Then <p is called a bump function. 

Definition 5.2 (Class Gp,<p)- Fix p G [1, oo) and let ip : R — ► R be a 
given bump function. Let .M be the set of all signed measures p, on 
(R, £>(R)) such that \p\ is a-finite. Define 

M PtV = {p e M : V G L[ oc (R, \p\) n L 1+ i(R, |/i|)}. 
Then for any p G M. v ^ define a function related to p by 

I(o,x] l P d ^^ forx>0, 
J(x,o]^ d ^ fori<0 > 
where (0, 0] = 0. Also define a set of jump functions 



A PiV = { A A (x) = ^ ^w{ai)X{ai}{x) ■ A = (a;)^ C R countable , 
i=i 

oo 

a» 7^ % if i ^ j, (Ai)~! C R and ^ | Ai|y?(ai) 1+ p < oo}. 

i=i 

Then we set 

G P ,<p = {g '■ g = c + g^ + for some c G R, 
p G «M Pj¥ , and A A G A p ^}, 

where the decomposition of g is unique, as we will see in Theorem 5.4. 
Moreover, denote the (p, ip) -variation of g G Q Pt(p by 

„ oo 

V p><p (g)= / <p 1+ pd\n\ + V |Ai|v?(aj) 1+ i 



14 
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Remark 5.3. The definition implies that any function g^ G Q v ^ is right- 
continuous and g M (0) = 0. To relax these restrictions, we add to the 
function g^ a constant c and a function A^, which can be used to 
alter the left- or right-continuity of g^ at the points of discontinuity or 
to add point discontinuities anywhere. For example, we can make g M 
left- continuous by choosing 



where A is the set of points of discontinuity of g^ . We see that g can 
have only a countable number of jumps, because otherwise /i({x}) ^ 
for uncountably many x G R, which is a contradiction to the a- 
finiteness of the measure Moreover, both g and g^ may jump at 
zero; the jump height of is then given by Lp(0)fi({0}), and the jump 
of g depends on (/?(0)/i({0}) and Aa(0). 

Theorem 5.4. The decomposition g = c + g M + A^ for functions g G 
Q p ,if is unique. 

Proof Take gi,g 2 £ G P ,<p such that = Q + g Mi + A^., i G {1,2}, 
and suppose that gi = g 2 - Now Ai U A 2 is countable and A^ = in 
(Ai U A 2 ) c . Let us take a sequence (xj) C (A\ U A 2 ) c such that Xj \ 
as j — > 00. Since g Mi is right-continuous and g^(O) = 0, we get that 
gi{xj) = Ci + g^ x (xj) — > Cj, and thus ci = c 2 . This implies that for 
x G (Ai U A 2 ) c we have g^(xQ) = g^ 2 (xo). Now let i G Ai U A 2 . 
Again we choose a sequence (xj) C (A 1 U A 2 ) c such that Xj \ rr as 
j — > 00, and by right-continuity of g Mi we get that g m (rr ) = g^ 2 (xo). 
Thus g w = (7^ 2 everywhere, and also A Al = A A2 . □ 

Theorem 5.5. Functions of bounded variation are a special case of 
functions in Q p ^, i.e. BV C G P ,<p- 

Proof. Let g G NBV and let fi BV be the signed measure related to 
g. At the points of continuity of g we have g = g(0+) + g^ 1 , where 
the measure fi is chosen such that dfi = dfi BV /ip on (0, 00) and dfi = 
—dfj, BV /ip on (—00, 0]. Here is a-finite by the finiteness of \i BV and 
the properties of ip. It also holds that 



Now let g G BV . Then g = g + c for some g G NBV and cG Rat the 
points of continuity of g, thus satisfying g = g(0+) + c + g^ . At the 
points of discontinuity we correct this by choosing A^ such that A is 
the set of the points where g is not right-continuous, and the values A 
correspond to the jump heights of g. Then 




g^{x—) — g^{x) for x G A, 
elsewhere, 



VWl 1+1 (R,M) 



< V(g)p+ 1 < 00. 



■oo 
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5.1. General Approximation. As before, let X and X be random 
variables. We define a function tp that connects the random variables 
with their tail behavior. 

Definition 5.6. Take two strictly positive monotone functions 

<p+: (0,oo)-(0,l] 

and 

p": (-oo,0] -(0,1] 
for < 9 < 1 with properties 

<Pq{K) — as if — oo, 
(Pff(K) — as if — — oo, 

P(X > if) V P(X > if )1 ° < ^(if) for if > 

and 

P(X < if) V P(X < if)] < <pj(K) for if < 0. 
Then we define a bump function <p X ' X : R — (0, 1] by 

' -g(K) ifif<0. 



Now the main result is the following convergence theorem for func- 
tions in the class Q Ptip associated with the function (pf' X . 

Theorem 5.7. Let X and X be random variables such that X has a 
bounded density. Suppose that < 9 < 1 and let ip e ' be a function 
as in Definition 5.6. If 1 < p < oo and g G Q p ^ x ,x , then for all 

1 < q < oo we have 



E\g(X) - g(X)\ p < 3 • 2 P (sup/x)^ (V p ^ x {gf 



X -X 



9(1-9) 
9+1 



Proof. Let g e G x,x- Then by definition g = c + g^ + and 

g(X)-g(X) < g^(X)-g^X) + A A (X) - A A (X) 

p v 

Now we can compute 



x,x. 



(0,x] 



z)dfx(z) X[o,oo)(x) 



/ X[z,oo)(x)cpf' x (z)dfi(z) 

J (0,oo ) 
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and similarly 



(x,0] 
f 

(-oo,0] 



V f' x {z)dfi{z) x(-oo,o)(x) 



X(-oo,z)(x) V f' x {z) dfi(z). 



Thus 



g»{X)-g»{X) 
+ 



< 



+ 



< 



+ 



<f(X) Xlo ,oo)(X) - g»(X) X[0 ,oo)(X) 
g»(X)x { -oo,o)(X) -^(^)X(-oo,o)W 

/ \X[z,oo)(X) -X[z,oo){X)\^f' X {z)d\fl\{ 

J (0,00) 



(-00,0] 



(0,00 

f 

(-00,0] 



X[ 2 ,oo)(^) -X[,,oo)(^) y fl ' (z) d\fi\(z) 

V 



X(-oo,z)(X) ~ X(-oo,z)(X) tpf X {z) d\n\(z). 

V 



Denote by ip(X, X) the error function from Theorem 3.4, i.e. 

X -X 



1 

9+1 



i{>(X,X) :=3(sup/x)^ 

and notice that a A b < a l ~ 6 b e for any a, b > and < < 1. Since 

E|x[x,oo)(X) - x^oo)(X)| = P(X > K,X< K) + P(X <K,X>K) 

< 2(P(X > K) V P(X > X)), 

it follows from Theorem 3.4 that, for K > 0, 

E|x[x,oo)PO - X [^oo)W| < A 2 [P(X > X) V P(X > K) 

< ^(X,X) l - 6 2 6 \P(X > K)VP(X > K) 
<2^(X,Xy- e cpt(K), 

where < 9 < 1. In a similar way we get for K < that 

nX[K,oo){X) - X [K,oo)(X)\ < 2 6 iP(X,X) 1 - e tpe(K), 
so we can write for K G R that 

nX\K t00 )(X) - X [K,oo){X)\ < 2^(X : Xy- e <pf' x (K). 
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This gives an estimate for X[z,oo)(X) — X[z,oo)(X) , and the same es- 



timate holds for 



X{-oo,z)(X) - X(-oo,z)(X) 



by the observation in Re- 



mark 3.5. Therefore 



p Jr 



e / l-e . , 9(1-9) 

<2p / 3 v (sup fx) (q+1)p 
Jr 



X -X 



(9 + l)P X,X( 71 \, \ 

ip e ' (z) pd\fi\(z) 



<3i(sup/ x )I^ / ip x ' x {z) 1+ l>d\ii\{z) 
Jr. 



X -X 



(?+l)p 



It remains to show a similar estimate for the jump function A A . This 
can be done by the same argument as in the case of bounded variation, 
namely by writing 



where 



Aa(x) = / X(-oo,z](x) -X(-oo,z)(x)dv(z), 



R 



v = ^2\ iV f' x {ai)5 at 



i=i 



and S a is the Dirac measure in a. Then by arguments similar to the 
first part of the proof and Remark 3.5 we get 



A A (X)-A A (X) 



< 


Ir 


p 





\X(-oo,z](X) - X(-oo,z](X)\d\is\(z) 



R 



< 



+ 



R 



\X(-oo,z)( X ) - X(~oo,z)(X)\ d\u\( Z ) 

X(-oo,z](X) - X(-oo,z](X) d\is\(z) 



X(-oo,z)(X) - X(-oo,z)(X) d\u\(z) 

V 



9(1-9) 



.1=1 



< 2 -3p(sup/ x )t5+iF / cpf' x (z)p d\is\(z) 

J R 

/ oo 

<2-3?(sup/ x )f™ [Y^WWe^iaiY^ 
so finally we get 



g(X)-g(X) <2.&(supf x )$&V x ,i(9) 

p P,<P e 



X -X 
X -X 



9(1-9) 
(9 + l)p 



9(1-9) 
(9 + l)p 



X -X 



9(1-9) 
(9 + l)p 



□ 
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6. Analysis of the Class Q p ^ 

We study the class Q P)lf> with the underlying function tp = tpf ,X - 
This function depends on the approximation X, and our first task is 
to handle this dependence. We show in Lemma 6.2 that we can choose 
the function ip e ' such that it decays faster than any polynomial, and 
then we prove in Theorem 6.3 that with this choice, the class Q X x 

contains all polynomials. Then we apply the results to solutions of 
SDEs, and collect our knowledge in the main result, Corollary 6.4. 

Lemma 6.1. Suppose that if and ip are bump functions. 

(i) If<p<ijj, then Q Pii> C Q Pjip . 

(ii) If g G Q p ^ and <p <ip, then V PtV> (g) < V Pti ,{g). 

Proof. First we show (i). Let g G Q p ^ and /i^ be related to g, i.e. 
g = c + gV* + A^. We choose a measure fj, v such that d/i^ = (ip/ip)d/i^, 
which implies that d\n v \ = / ip)d\n$\ and \fi v \ is a-finite. Now we get 
for x > that 

g^(x)= i/j(z) dm(z) = <p(z) dn v (z) 

J{0,x] J(0,x] 

and similarly for x < 0. The integrability conditions are satisfied, since 
ip(z) d\ii, v \(z) = / ip(z) d\^\{z) < oo 

A J A 

for all A CC R and 

/ (p(z) 1+ p dln^Kz) = / ijj(z)ip(z) rd\fi t I (z) 
Jr. Jn 

< / %l)(z) 1+ p d\n^\{z) < oo. 

J R 

The representation of the jump part A^ changes correspondingly in 
the change of measure, i.e. we set 



Xf — X 



<p(Oi) 

and see that 

oo oo 

^|Ari^) 1+ l<^|Af|^) 1+ i<oo. 

i=l i=l 

This proves the assertion (i), and (ii) follows by a similar argument. □ 

Lemma 6.2. Suppose that X and X are random variables such that 
X G rip6[i,oo) L P , and suppose there exists C = (C p ) pe [i )00 ) C (0, oo) 



such that 



X -X 



< C p for all p G [1, oo). let 9 G (0, 1). Then we 

p 

X ' X n„.„U +h„+,r, X ' X — ,„X 



can choose the function ip e ' such that <p e ' = Pe g, where the function 
(Pqq is a bump function that decays faster than any polynomial. 
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Proof. The triangle inequality gives that X G L p and 

+ m P <c P +\\x\\ p . 

Thus by Chebychev's inequality we have for all A > that 

E|X| p 
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X 


< 


X -X 




p 





P(\X\ > A) < 



and 



P(|^| > A) < 



e\x\? k (c p + \\x\\ p y 



\p - \p 

So we have a polynomial tail estimate for X and X that depends only 
on the constants C p of the L p -estimates, not directly on X. This implies 
that 



P{\X\ > A) V P{\X\ > A) 



P e[i,oo) \ e P 



< inf 

e[i,c 



for A > 0. For A < we define ip^(\) := ^c'ed^D' an d Pc'di®) := 

X 

1. The function (p c ' d satisfies the monotonicity properties of a bump 
function, but is not necessarily strictly positive. However, if we take a 
bump function ip and define 

y X X X 

then ifcj is a bump function suitable for the choice of (f d ' . Since 
clearly decays faster than any polynomial and we can choose if)(\) = 
e~l A l, we see that (p^ e also decays faster than any polynomial. □ 

Let V be the set of all polynomials from R to R. Then we have the 
following: 

Theorem 6.3. Suppose that (p is a bump function that decays faster 
than any polynomial. Then V C Q p ^ for all p G [1, oo). 

Proof. Let p G [l,oo) and suppose that g G V. Then by the funda- 
mental theorem of calculus we have for x > that 



A 1 =: y^(A) 

Xfi, 



g(x) = g(0)+ I g\z)dz 
Jo 



and for x < that 



g(x) = g(0)- I g\z)dz. 

J X 

Thus by defining c = g(0) and a signed measure \i such that 



= sgn(z) 



we have that is a-finite and the representation g = c + g^ holds. 
Now g' also has only polynomial growth, say |<7'(x)| < C(l + |:r| s ) for 
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s > 1. But tp decays faster than any polynomial, so we have (p(x) < 
C\ x \-p('+ 2 ) A 1 and 

V p,M = Jr</ + MH = J n (pH z )W(z)\dz 

< CC~v / R (|^r {s+2) A 1)(1 + \z\ s ) dz < oo, 

which implies that g G G P ,ip- D 

Let us now come back to the SDE (2.1) and summarize our knowl- 
edge: 

Corollary 6.4. Let p G [1, oo). Suppose that X is the solution of the 
equation (2.1), X T has a bounded density, and X^ is an approximation 
of X T such that 

\\X T - Xt\ p < C p |7r| 7 
for some constants 7 > and C p > 0. Then for any < e < 7 we 
have for 6 = V^ce according to Lemma 6.2 and g G Q pip x T that 



E 



g{X T )-g{X®\* < 3 - ^(sup/^) 1 -^ (V ^{g))* &~j M 



7-e 



Especially, V C Q x T . 

Proof. By Lemma A. 2 we have that X? G C\ p e[i 00) A" so by Lemma 

Xt X 71 x X 

6.2 we can choose ip e ' T = y? C g, wnere (f c T e is a bump function with 
decay faster than any polynomial. Now using Theorem 5.7 we get for 
any q G [1, 00) and 9 G (0, 1) that 

E\g(X T ) - g(X*)\ p 

(2(1-0) / \P li 1 - 9 ) 79(1-9) 

<3-2^(su P / Xt )— (V ptV xr{g)) C q ^ |tt|— . 

Let < e < 7. Now choose q = ^ — 1 and let = 1/q. Note that 
g > 1 since e < 7. Then 

g(l - 0) _ g - 1 _ 1 £ 
g + 1 g + 1 7 
and thus we get for all g G Q x T that 



E 



7-e 



Moreover, by Theorem 6.3 we have that V C Q x T . □ 

Remark 6.5. In Corollary 6.4 the function (p e ' T depends on the dis- 
tribution of X£ and is replaced by the uniform bound ipQ T e . However, 
when considering convergence rate we are looking at partitions with 
small mesh size. Thus if approximating random variables X£ corre- 
sponding to partitions with large mesh size had heavy tailed distribu- 
tions, the use of the uniform bound could unnecessarily narrow down 
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the class of functions. Therefore in such 8b C£lS6 it would be better to 
take more delicate approach and study the result 



7-e 



E\g(X T ) -g{X%)\* < 3 • 2"(sup f XT f~i (v * T ,* f (g)Y c\,? |tt| 

Corollary 6.4 now gives convergence rates for both Euler and Milstein 
schemes: 

Corollary 6.6. Let 1 < p < oo. Then for < e < \, 9 — and 

g G Q x T that 

E\g(X T ) - g(X*)\> < 3 • 2^(sup/ XT ) 1 - 2£ (V ^ (g))" C{j* . 
and similarly for < e < 1 and 9 = ^ we have that 

E\g(X T ) - g{X?)\* < 3 • 2^(sup/x T ) 1 - £ (V ^{g))* C\rf ^ ■ 
Especially, the statements hold for any g G V . 

Example 6.7. Let us generate a jump function by choosing the measure 
ji to be a sum of Dirac measures, 

m X T( n \ fe ' 

where a^, a& G R for all k G Z and 7^ a; for k ^ I. Then from the 
integrability condition for /x we see that g <E G x T if 

(6.1) |« fc |v?c,eK) P < °°- 

Therefore the result of Corollary 6.4 holds for jump functions with 
jumps controlled by the decay of the function ip^ T g in a way that the 
condition (6.1) is satisfied. 

6.1. Euler scheme. In the case of the Euler scheme we can again 
use our knowledge about constants to get more explicit results for the 

decay of the function ip e ' T . Let us recall the following result from 
Bouleau and Lepingle [1]: 

Lemma 6.8 ([1, Ch. 5, Lemma B.1.2.]). For 1 < p < 00, 

sup\X?\ G LP 

t<T 

and there exist M(xq,T,Ct) > such that 



sup \X } 

t<T 



<e Mp \ 



_ e 
v 
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Proof. By Theorem A.l and Lemma A. 2 in the Appendix we get 



sup \Xf\ 



t<T 



< 


sup|X t -Xf| 


+ 


sup \X t - x \ 


+ \%o\ 




t<T 


P 


t<T 


p 



<(2VT+\x \)e Mp2 , 
and we absorb the constant (2VT~ + \x \) into the constant M. 



□ 



Xt,XS 



Theorem 6.9. We can choose the function tp e ' T in a way that 
^ t ' x t < (pfs, where ipf^ is a bump function such that 

(i) if the functions a and b are bounded, i.e. \a\, \b\ < M, we have 



, 2aPt if z > max(i + MT, 0), 

_e(z-(x -MT)f 
2M 2 T 



if z < min(a;o — MT, 0), 
elsewhere, 



(ii) if the functions a and b are Lipschitz, then there exists z > 1 
such that we have 

[1 if \z\ < z , 

where M = M(x , T, C T ) > 0. 

Proof, (i) We consider the Euler approximation with n time nodes in 
the integral form (2.2). If we denote 

n-l 

L u ■= ^ff(4,^)X[t t ,t H1 )(tt), 

k=0 

then by the boundedness of a and the Novikov condition 



(lu 



is a martingale for any a > 0, and EM ( = 1. Thus by Chebychev's 
inequality we have for A > 1 that 

P (e a $o ' L u dW u -2§- / T 'l£du > < _J_ 

By taking logarithm this implies 

(pT 2 pT \ 

a J L u dW u -^-J L 2 u du>\j<e- X 



for A > 0. Since 



L 2 u du < M 2 T, 



we get 



_ , . T „ TT A aM 2 T 



<e'\ 
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which we can reparametrize to get 

pQT \ u dW u >X^<e^^ 
for A > aM 2 T/2. Now we can choose a = A/(M 2 T) to get 

P Qf L u dW u > A^ < 
for A > 0. A similar proof with L u = —L u shows that 
P ( / /.„ d\Y., < x] < e~^r 



o 



for A < 0. Therefore, for A > max(x + MT, 0) 

P (X? > A) < P + jT L u dW u + MT>\^j 



(\-(x +MT)) 2 



< e 2M 2 T 

and for A < min(xo — MT, 0) 



P (X| < A) < P + f o L " dW ^ -MT <\^j 

(\-(x Q -MT)) 2 
<C e 2M 2 T . 

Obviously a similar proof works for the random variable X T instead 

of X-f , so by the definition of (p e ' T the assertion follows. Moreover, 
to get a bump function we choose the upper bound to be one on the 
interval [min(x — MT, 0), max(x + MT, 0)]. 

(ii) If a and b are Lipschitz, then we know from Lemma 6.8 that 

\\xm <e Mp \ 

II 1 lip — 

where the constant M > depends on x , T and Ct- Now by Cheby- 
chev's inequality we have for A > that 

P(|*#| > A) < -LfL < 



Choose 3Mp 2 = log A for A > A = e 3M . This gives 

_ (log A) 1/2 
P (3M)V2 ' 

and thus for A > Ao we get 

P(|X#| > A) = ^ = A-i' = A-i^) 1/2 . 

VI Tl- I Xp 

Again the same proof works for the term P(|X T | > A) because of 
Lemma A. 2 in the Appendix. □ 
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Theorem 6.10. Let c > 0. // the functions a and b are bounded, and 

g{z) = e c ^ z/ < 7 < 2, 

or 

g(z) = e c|z| with c < 9/p, 

then g e Q x T . 

Proof. Since g is not differentiable at zero, define g(x) := g'(x) if x ^ 
and g(0) := 0. By choosing a signed measure 

sgn(z)#(z) 

we get that is a-finite and the representation g = g(0) + g^ holds. 
The definition of the class Q x T gives the condition 

/ {^ E {z)Y\~g{z)\dz= [ {^(z^e^c^z^dzKoc, 

which is by Theorem 6.9 satisfied, because the singularity at zero for 
< 7 < 1 is not too strong, and integrability is determined by the 
parameters 7, c, 9 and p as proposed in the formulation of this Theorem. 
Similarly we see that the local integrability condition is satisfied. □ 

7. Lower bound 

In this section we find a solution X\ (i.e. T = 1) of an SDE of the 
type (2.1) such that it gives a lower bound for the approximation rate of 
the Euler scheme in Theorem 3.7. This is achieved by choosing X t = St, 
the geometric Brownian motion. Let St = e Wt ~ t l 2 for t G [0, 1], so that 
S is a solution of 

S t = l+ f S s dW s 





and let U n := S E — S, where S E is the Euler scheme as defined in (2.2) 
corresponding to the equidistant partition of [0, 1], i.e. n = (i/n)™ =0 . 

Lemma 7.1. We have (W, y/nll n ) =^ (W, U) in the Skorohod topol- 
ogy, where U is the strong solution of the equation 

(7.1) U t = ! U s dW s -^= [ S s dB s 



V2, 

and B is a standard Brownian motion independent of W . 

Proof. The statement is an immediate consequence of a result by Jacod 
and Protter, [10, Corollary 5.4]. □ 

Theorem 7.2. There exists K > such that 

liminf sup E|x[*r,oo)(<Si) - X[K,ao)(Sf)\ > 0, 
n^oo K y Ko 

where Sf is the equidistant Euler approximation of Si. 
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Remark 7.3. Theorem 7.2 states that the convergence rate ^ ~ e f° r the 
Euler scheme obtained in Theorem 3.7 and consequently in Theorem 
4.5 and Corollary 6.6 is optimal up to the factor e, i.e. any rate 7 > | 
leads to a contradiction with the statement of Theorem 7.2. 

Proof of Theorem 7.2. Let us consider the setting of Lemma 7.1 and 
the process U defined by the equation (7.1). If Lq — a.s., then for all 
t G [0, 1] we have Ut — a.s., which leads to a contradiction. Therefore 
P (U 1 > 0) > or P (U 1 < 0) > 0. If P (E/i > 0) > 0, then there exist 
e G (0, 1], 5 > and K > 1 + K Q with K > such that 

P(Si G [fT-l,fO, [q >e) =<?. 

The case P (Lq < 0) > can be treated in a similar way by chang- 
ing the condition Lq > e to Lq < — e. By Lemma 7.1 we know 
that (W 7 ", s/nU n ) =>• (W, U) in the Skorohod topology. This implies 
that (Wi,y/nU™) =>■ (W^Lq), since the projection mapping 7iq, i.e. 
the mapping a \— > a(l) for a process a, is continuous in the Sko- 
rohod topology. Because the function is continuous, we have 
(Si,y/nU?) (Si,l/i). Therefore 

liminf P (5i G [if - 1, if), >/n[£f - Si] > e) 
= liminf P K G (X- VnC/r > e) 

n— »oo 

> P (Si G (K-1,K),U! > e) 
— P (Si G [K — 1, X), Lq > e), 
and we see that there exists n > 1 such that for all n > no 

pUe[K-l,K), [s*- Sl ]>-Z=)> 5 -. 



Assume a partition K — 1 = K™ < K™ <••• < K™ = K. Then 

5 



sup P (s l G p-!, AT), [5f - S l ] > -^=) 

l=l,...,m V V n J 



> 

2m 



Now choose the partition (iff 1 )™! to be equidistant with 
(7.2) - < -i=. 



Then there exists /o G {1, . . . , m} such that 
5 



< 

2m 



<P{S 1 < K%, 5f > . 



Let m = [-^/n/^], which satisfies the condition (7.2) for the mesh size. 
Hence 

£ V( Sl <K Z Sf>KO 
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< E|x[X-,oo)(5'l) - X[K%,oo)(Si)\. 

Since \y/n/e] < we have 



E|X[Jfi",oo)('S , l) - X[K™,oo){Sf)\ > 



> 



5e 



Therefore 



Se 



Vn sup E|x[^,oo)(^i) -X[x,oo)(^i )| > -j 



K>K 



for all n > n , which implies the assertion. 

Appendix A. 



□ 



Here we prove the following Theorem from the book of Bouleau and 
Lepingle, [1, pp. 275-276]. The proof is given in the book, but without 
computing the constant explicitly. 

Theorem A.l. If the assumptions (i)-(iii) in section 2 hold, and 1 < 
p < oo, then 

sup \X t -X t E \ <e Mp2 \n\^, 

0<t<T 

where the constant M > depends at most on x , T and Ct- 

For the proof we need the following Lemma: 
Lemma A. 2. For 1 < p < oo and < s < t < T , 



sup \X U - X s 



s e 



Mp 2 



where M > depends at most on x , T and Ct- 

Proof. Without loss of generality we can suppose that p > 2. For fixed 
s e [0, T] define 



F(t) 



sup \X U - X s 

u£[s,t] 



for all t G [s,T]. Then by the Burkholder-Davis-Gundy inequality (as 
in [1, p. 269]) and the linear growth condition we get 



F(t) < 2 
< Sp 



a(u,X u ) dW u 



\a(u,X u )\ 2 du 



< C T (8p+Vt^) 



- I \b(u,X u )\du 

+ Vt - s 
(l + \X u \) 2 du 



\b(u,X u )\ 2 du 



CONVERGENCE RATES FOR SDES 



27 



< c 



< c 



< c 



(1 + \X s \fdu 



Vt=s\l + \X s \\\ p + 



Vt^~s(l + \\X s \\ p ) 



\X U -X s \) 2 du 



sup \X V — X s \ 2 du 
s ve[s,u] 



sup \X V — X s 

vE[s,u] 



du 



= Cy/t=l(l + \\X S \\ ) +C[ / F{u) 2 du) , 



where C = Cx(Sp + \/t — s). Thus we have 

F{tf < 2C 2 (S P +Vt^~s) 2 (t-s)(l + \\X s \\ p ) 2 

+2C^(8p + Vt^sf F{ufdu 

and by the Gronwall lemma we get 

F(t) 2 < 2C 2 (8p+Vt^) 2 (t - s)(l + lx s \\ p ) 2 e 2C T^ + ^ 2{t - s \ 
By taking the square root and choosing a suitable constant M we get 

F(t) < MpVt^~s(l + \\X s \\ p )e M v 2 , 
where M = M(T, Ct)- In particular, the above estimate gives that 

\\X s \\ p <\x \ + \\X S - X \\ p < \x Q \ + Mp^(l + \x \)e Mp \ 
so if we redefine the constant M, we get 

< \ft~- s ( 



sup \X U - X s 

ue[s,t] 



Mp 2 



where M = M(x ,T,C T ). 

Proof of Theorem A.l. Suppose that p > 2. Now define 



□ 



F(t) 



sup \X S - X 

s<t 



Here X E is the Euler scheme related to the equidistant partition it, 
and is defined for continuous time by formula (2.2). Then by the 
Burkholder-Davis-Gundy inequality ([1, p. 269]) we get 



F(t) < 8p 



r-t n-l > 

X[t k ,t k+1 )(u) du 
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+ Vt 



n-1 



E Hn,X u ) - b(t k ,X t E k )\ 2 X[tk , tk+l) (u) du 



k=0 



'n-1 



^l 6 ( u ' X 0-K*fc>^)| p X[t fc> t fc+ i)(«)j du 
By the conditions in Section 2 we have that 



\a(u,X u ) -cr(t k ,Xf k )\ < \a(u,X u ) -a(u,X t 
+ \a(u,X tk )-a(t k ,X tk )\ + \a(t k ,X tk )-a(t k ,Xfj\ 

< c T (\x u - x tk \ + (i + \x tk \)\u -t k r + \X tk - xf k \) , 

and a corresponding inequality holds for the function b. Thus by 
Lemma A. 2 we have for u > t k that 

\\a(u,X u )-a(t k ,Xf k )\\ p 

< C T (\\X U - X tk \\ p + 1(1 + \X tk \)\u - t k n\ p + \\x tk - x t E k \\ p ) 

< C T (e Mp2 \u - t fc |3 + (l + \x Q \ + Vhe Mp2 ) \u - t k \ a + F(u)) 

< C T ((1 + |x | + >/T)e Mp2 (|7r|5 + |7r| Q ) + F(«)) , 

and again a corresponding inequality holds for b. Denote C(x ,T) = 
1 + \x \ + \/T. Now we can continue our estimate for F(t) to get 

F(t) < C T {8p + Vt) • 

/ Y,\ c ^ T ) eMp2 (\^ + \^n+ F ( u ) x[t kM { 



2 \ 5 

du 



<C T (8p+VT) ( / (C(x ,T)e M ^(|# + M'\) + ./'(>/)) f/» 

< c T (8p + 

< e Mlp2 (|7r|5 + |tt| q ) + M lP ^jf F{uf 







i -i 

2 




/ F(ufdu 






Jo 
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where M 1 = Mi(xq,T,Ct) does not depend on p. Taking a square we 
get 

F{tf < 2e 2Mlp2 (|vr|5 + \n\ a ) 2 + 2M 2 p 2 f F(u) 2 du, 

Jo 

and thus Gronwall's Lemma gives 
F(t) 2 < 2e 2Mlp2 (|7r|^ + \n\ a ) 2 e 2 ^ T = 2e^ +2M " T)p2 (\ir\^ + \ir\ a ) 2 . 

By taking square root on both sides and recalling the assumption a > | 
this gives 

F(t) < V2e^ m+M ^ p2 (\7r\^ + \ir\ a ) < e M ^\^, 

where M 2 = M 2 (x ,T,C T )- The case 1 < p < 2 follows from the case 
p = 2 by redefining the constant M 2 . □ 

References 

[I] Nicolas Bouleau, Dominique Lepingle, Numerical Methods for Stochastic 
Processes. Wiley, 1994. 

[2] Vlad Bally, Denis Talay, The Law of the Euler Scheme for Stochastic Dif- 
ferential Equations: I. Convergence Rate of the Distribution Function. Probab. 
Theory Related Fields 104 (1996), no. 1, 43-60. 

[3] Vlad Bally, Denis Talay, The Law of the Euler Scheme for Stochastic Dif- 
ferential Equations: II. Convergence Rate of the Density. Monte Carlo Methods 
Appl. 2 (1996), no. 2, 93-128. 

[4] Colin Bennett, Robert Sharpley, Interpolation of Operators. Academic 
Press, 1988. 

[5] Maria Emilia Caballero, Begona Fernandez, David Nualart, Esti- 
mation of Densities and Applications. J. Theoret. Probab. 11 (1998), no. 3, 831 
- 851. 

[6] Avner Friedman, Partial Differential Equations of Parabolic Type. Prcnticc- 
Hall, 1964. 

[7] Norbert Hofmann, Thomas Muller-Gronbach, Klaus Ritter, The 
Optimal Discretization of Stochastic Differential Equations. J. Complexity 17 
(2001), No.l, 117- 153. 

[8] Norbert Hofmann, Thomas Muller-Gronbach, Klaus Ritter, Linear 
vs. Standard Information for Scalar Stochastic Differential Equations. J. Com- 
plexity 18 (2002), 394 - 414. 

[9] Norbert Hofmann, Thomas Muller-Gronbach, On the Global Error of 
ltd- Taylor Schemes for Strong Approximation of Scalar Stochastic Differential 
Equations. J. Complexity 20 (2004), 732 - 752. 

[10] Jean Jacod, Philip Protter, Asymptotic Error Distributions for the Euler 
Method for Stochastic Differential Equations. Ann. Prob. 26 (1998), no. 1, 267- 
307. 

[II] Ioannis Karatzas, Steven E. Shreve, Brownian Motion and Stochastic 
Calculus, Second Edition, Springer- Verlag, 1991. 

[12] Peter E. Kloeden, Eckhard Platen, Numerical Solutions of Stochastic 

Differential Equations. Springer- Verlag, 1992. 
[13] Thomas Muller-Gronbach, Strong Approximation of Systems of Stochastic 

Differential Equations. Habilitation thesis. Darmstadt, 2002. 



30 



R. AVIKAINEN 



[14] Thomas Muller-Gronbach, The Optimal Uniform Approximation of Sys- 
tems of Stochastic Differential Equations. Ann. Appl. Probab. 12 (2002), no. 2, 
664 - 690. 

[15] Thomas Muller-Gronbach, Optimal Pointwise Approximation of SDEs 
Based on Brownian Motion at Discrete Points. Ann. Appl. Probab. 14 (2004), 
no. 4, 1605-1642. 

[16] Kyoung-Sook Moon, Anders Szepessy, Raul Tempone, Georgios E. 
Zouraris, Convergence Rates for Adaptive Weak Approximation of Stochastic 
Differential Equations. Stoch. Anal. Appl. 23 (2005), no. 3, 511-558. 

[17] David Nualart, The Malliavin Calculus and Related Topics. Springer- Verlag, 
1995. 

[18] Walter Rudin, Real and Complex Analysis, Second Edition. McGraw-Hill, 
1966, 1974. 

Department of Mathematics and Statistics, P.O. Box 35 (MaD), FIN- 
40014 University of Jyvaskyla, Finland 
E-mail address: avikain@maths.jyu.fi 



